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Abstract 

The present status of our ongoing systematic study of the discovery 
potential of QCD-instanton induced events in deep-inelastic scatter¬ 
ing at HERA is briefly reviewed. We emphasize our recent progress 
in predicting the cross-sections of instanton-induced processes. Our 
finalized predictions include a dramatic improvement of the residual 
renormalization-scale dependencies and the specification of a “fidu¬ 
cial” kinematical region in the relevant Bjorken variables extracted 
from recent lattice simulations. Published upper limits on instanton- 
induced cross-sections based on single observables in the final state 
such as the flow of strange particles and the multiplicity distribution 
of charged particles are already of the order of our estimate. Thus, a 
decisive search for instanton-induced events in deep-inelastic scatter¬ 
ing at HERA, based on a multi-observable analysis, seems feasible. 


*Talk presented at Quarks ‘98, 10th International Seminar on High Energy Physics, 
Suzdal, Russia, May 1998; to be published in the Proceedings. 



1 Introduction 


In this contribution, we briefly review the present status of our ongoing 
systematic study 0-^] of the discovery potential of QCD-instanton induced 
events in deep-inelastic scattering (DIS) at HERA. 

Hard scattering processes in elementary particle physics are successfully 
described by the Standard Model of strong (QCD) and electro-weak (QFD) 
interactions in its perturbative formulation. However, perturbation theory 
does not exhaust all possible hard scattering processes. 

Instantons [0 , fluctuations of (non-abelian) gauge fields representing topol¬ 
ogy changing tunnelling transitions, induce interactions which are absent 
in conventional perturbation theory. In accord 0 with the general chiral 
anomaly relation, instantons give rise to (hard) processes in which certain 
fermionic quantum numbers are violated, notably, chirality (Q 5 ) in (massless) 
QCD and baryon plus lepton number (B -j- L) in QFD. 

Implications of QCD-instantons for long-distance phenomena have been 
intensively studied in the past, mainly in the context of the phenomenological 
instanton liquid model P and of lattice simulations |T^. Yet, a direct exper¬ 


imental verification of their existence is lacking up to now. An experimen¬ 
tal discovery of such a novel, non-perturbative manifestation of non-abelian 
gauge theories would clearly be of basic significance. 

The deep-inelastic regime is distinguished by the fact that here hard 
QCD-instanton induced processes may both be calculated |TI|flJl[| within in- 
stanton-perturbation theory and possibly detected experimentally As 

a key feature it has recently been shown [Q], that in DIS the generic hard 
scale Q cuts off instantons with large size p S> over which one has no 
control within perturbation theory in the instanton background. 


2 Cross-Sections 

The leading instanton (/)-induced process in the DIS regime of e^P scatter¬ 
ing is displayed in Fig. |l[ The dashed box emphasizes the so-called instanton- 
subprocess with its own Bjorken variables, 

Q'2 ^ > Q. ^ 

2 ,p ■ q 

The /-induced cross-section in unpolarized deep-inelastic e^P scattering. 
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Figure 1 : The leading instanton-induced process in the DIS regime of e^P 
scattering, violating chirality by AQ5 = 2 nf. 


daep, can be shown to factorize [|^ (in the Bjorken limit, = —q^ large) 
into a sum of differential parton-parton luminosities, dCp,p, and J-subprocess 
cross-sections, a pip, 


da, 
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( 2 ) 


Here p' = q',q' denotes the virtual quarks entering the J-subprocess from 
the photon side and p = q,q,g denotes the target partons. The differential 
luminosity dCp,p, accounting for the number of p’p collisions per eP collision, 
has a convolution-like structure |^, involving integrations over the target- 
parton density, fp, the 7*-flux, P^*, and the known 0 flux PpP of the parton 
p' in the /-background. 

In Eq. (]^), the J-subprocess total cross-section apip contains the essential 
instanton dynamics. We have evaluated the latter by means of the opti¬ 
cal theorem and the so-called JJ-valley approximation for the relevant 
q'g q'g forward elastic scattering amplitude in the II background. This 
method resums the exponentiating final state gluons in form of the known 
valley action and reproduces standard /-perturbation theory at larger 
II separation a/^. While a semi-quantitative evaluation of a pip may also 
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be performed by means of standard J-pertnrbation theory, the nse of the II- 
valley method allows to extend onr calcnlations to a somewhat larger range 
in x'. 

Corresponding to the symmetries of the theory, the instanton calcnlns in¬ 
troduces at the classical level certain (undetermined) “collective coordinates” 
like the / (/)-size parameters p (p) and the II distance B?/pp (in units of 

the size). Observables like must be independent thereof and thus involve 
integrations over all collective coordinates. Hence, we have generically, 

OO C50 

^P'l = J J dpD{p) J d^R ... (3) 

0 0 

X Q-(P+P)Q' _ 


The hrst important quantity of interest, entering Eq. (|^), is the /-density, 
D{p) (tunnelling amplitude). It has been worked out a long time ago |8|JT3[ 
in the framework of /-perturbation theory: (renormalization scale pr) 


D{p) = d 


^ — /3o + 



exp 


27r Apprf 


(h'j’ 

(A - 12/do), 




( 4 ) 


( 5 ) 


in terms of the QCD /d-function coefficients, /do = 11 — |n/, /di = 102 — 


38 


Uf. In this form it satishes renormalization-group invariance at the two- 


loop level [0. Note that the large, positive power 6 of p in the /-density 
would make the integrations over the /(/)-sizes in Eq. (^) infrared divergent 
without the crucial exponential cut-off 0 e-(p+p)Q arising from the virtual 
quark entering the /-subprocess from the photon side. 

The second important quantity of interest, entering Eq. (|^), is the II- 
interaction, — 1. In the valley approximation, the //-valley action, 

S(/D = 

is restricted by conformal invariance to depend only on 
the “conformal separation”, ^ /pp + p/p + p/p, and its functional form 

is explicitly known 0- It is important to note that the interaction between 
I and I remains attractive for all separations A specihcally, the //-valley 
action decreases monotonically from 1 at inhnite conformal separation to 0 
at = 2, which corresponds to /?^ = 0 and p = p. 

The collective coordinate integration in the cross-section (^ can be per¬ 
formed via saddle-point techniques. One hnds R* = (p*a/^* — 2,0) and 
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Figure 2: /-subprocess cross-section |^. 


p* = p*, where the saddle-point solutions p* and behave qualitatively 
as 

dvr 


* 

p ~ 


asQ'' 


e-2 = 


R* 

p* 



( 6 ) 


gluon, aljPg 


Thus, the virtuality Q' controls the J(/)-size: As one might have expected 
intuitively, highly virtual partons probe only small instantons. The Bjorken- 
variable x', on the other hand, controls the conformal separation between I 
and I: for decreasing x', the conformal separation decreases. 

Our quantitative results |]^ on the dominating cross-section for a target 
are shown in detail in Fig. both as functions of Q'^ (left) and 
of x' (right). The dotted curves in Fig. |, indicating lines of constant p* (left) 
and of constant R*/p* (right), nicely illustrate the qualitative relations (^ 
and their consequences: the Q' dependence essentially maps the /-density, 
whereas the x' dependence mainly maps the //-interaction. 

Compared to our earlier estimates based on the one-loop renormalization- 
group invariant form |@] of the /-density D{p), the residual dependence on 
the renormalization scale has now dramatically weakened (c.f. Fig. (left)). 
This important stabilization of our predictions originates from the use of the 
two-loop renormalization-group invariant improvement (^ ^) of D{p) from 
Ref. 0. 

Intuitivelely one may expect 00 Pr ~ l/{p) ~ Q'/Po = 0(0.1) Q'. 
Indeed, this guess turns out to match quite well our actual choice of the 
“best” scale, pr = 0.15 Q', determined by da^Pg/dpr ^ 0 (c.f. Fig. || (left)). 

Important information about the range of validity of /-perturbation for 
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Figure 3: Left: Renormalization-scale dependence of the J-subprocess cross- 
section for a target gluon []^. Right: Support for the validity of I- 
perturbation theory for the /-density D{p) from recent lattice data m. 


the /-density and the //-interaction, in terms of the instanton collective 
coordinates (p < Pmax,/?/p > (/?/p)min), can be obtained from recent (non- 
perturbative) lattice simulations of QCD and translated via the saddle-point 
relations (|^) into a “fiducial” kinematical region {Q' > |l[. 

In fact, from a comparison of the (one-loop) perturbative expression of the 
/-density (|[) with recent lattice “data” [|^ one infers [|^ semi-classical /- 
perturbation theory to be valid for p<pmax — 0.3 fm (c.f. Fig. I (right)). 
Similarly, it can be argued that the attractive, semi-classical valley result 
for the //-interaction is supported by lattice simulations down to a minimum 
conformal separation / min — 3, i.e. (/?*/p*)min — 1- The corresponding 
“fiducial” kinematical region for our cross-section predictions in DIS is then 
obtained as [Q 

(f < 0.3 fm; 1 f 0 ' > ~ 8 GeV; 

f > 1 P R' > XX 0.35, 

Fig. § displays the finalized /-induced cross-section at HERA, as func¬ 
tion of the cuts and obtained with the new release “QCDINS 

1.6.0” 0 of our /-event generator. For the following “standard cuts”, 

Cstd = x' > 0.35, Q' > 8 GeV, xbj > 10■^ 0.1 < pbj < 0.9, (8) 

including the minimal cuts (|^ extracted from lattice simulations, we specif¬ 
ically obtain 

'^HERA(^std) = 126 pb, (9) 
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Figure 4: /-induced cross-section at HERA []^. 


where the uncertainties result mainly from the experimental uncertainty in 
the QCD scale A. Hence, with the total luminosity accumulated by experi¬ 
ments at HERA, L = 0(80) pb“^, there should be already 0(10'^) /-induced 
events from the kinematical region (H) on tape. Note also that the cross- 
section quoted in Eq. corresponds to a fraction of /-induced to normal 
DIS (nDIS) events of 

= = ( 10 ) 

'^HERA (hstd) 

Thus, it appears to be a question of signature rather than a question of 
rate to discover /-induced scattering processes at HERA. Hence, we turn 
next to the final states of /-induced events in DIS. 

3 Searches at HERA 

A Monte-Carlo generator, QCDINS, for instanton-induced events in DIS, in¬ 
terfaced to HERWIG, has been developed |Q-0. The Monte-Carlo simulation 
essentially proceeds in three steps. First, quasi-free partons are produced by 
QCDINS with the distributions prescribed by the hard process matrix ele¬ 
ments. Next, these primary partons give rise to parton showers, as described 


7 




2.5 


> 

0 ) 

O 



2 

1.5 

1 

0.5 

0 



Figure 5: Lego plot of a typical /-induced event in the HERA lab system. 


by HERWIG. Finally, the showers are converted into hadrons, again within 
HERWIG. 

In Fig. ^ we display the lego plot of a typical /-induced event at HERA, 
as generated by QGDINS. Its characteristics directly reflect the essential fea¬ 
tures of the underlying /-subprocess and thus may be intuitively understood: 

After hadronization, the current quark in Fig. |l] gives rise to a current- 
quark jet. No further distinct qjets are expected, since the partons from 
the /-subprocess are emitted spherically symmetric in the p'p c.m. system 
(“/-c.m. system”). The gluon multiplicities are generated according to a 
Poisson distribution with mean multiplicity ~ l/ttg ~ 3. In view of 

the large required chirality violation AQ^ = 2n/ = 6, the total mean parton 
multiplicity is large, of the order of ten. After hadronization, we therefore 
expect from the /-subprocess a hnal state structure reminiscent of a decaying 
hreball: > 20 hadrons are produced, always including strange ones. They are 
concentrated in a “band” at hxed pseudorapidity p in the ( 77 , azimuth angle 
0)-plane. Due to the boost from the /-c.m. system to the HERA-lab system, 
the center of the band is shifted away from 77 = 0. Its width is of order 
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At] ~ 1.8, as typical for a spherically symmetric event. For x' ~ 0.35 and 
Q' ~ 8 GeV, the total invariant mass of the /-system, a/s' = Q'— 1, 
is expected to be in the 10 GeV range. All these expectations are clearly 
reproduced by our Monte-Garlo simulation. 

These features have been exploited by experimentalists at HERA to place 
hrst upper limits on the fraction of /-induced events to normal DIS events, 
in a similar kinematical region as our standard cuts (§): From the search of a 
excess in the “band” region, the HI Gollaboration could establish a limit 
of /[;2 = 6 %, while the search of an excess in charged multiplicity yields 
/lim = 2.7 % The limit from the charged multiplicity distribution has 

been further improved in Ref. to about 1 %. 

Since these experimental upper limits already range close to our esti¬ 
mate (11^), it becomes extremely interesting to investigate further possibil¬ 
ities to discriminate /-induced from nDIS hnal states. A dedicated multi¬ 
observable analysis is in progress, with the aim of producing an instanton- 
enriched data sample. Furthermore, strategies to reconstruct x' and Q'^ from 
the hnal state are being developped. 

Altogether, a decisive search for /-induced at HERA appears to be feasi¬ 
ble. 
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